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Abstract: Secondary electron cascades are responsible for significant ionizations in macroscopic sam- 
ples during irradiation with X-rays. A quantitative analysis of these cascades is needed, e.g. for assessing 
damage in optical components at X-ray free-electron lasers, and for understanding damage in samples exposed 
to the beam. Here we present results from Monte Carlo simulations, showing the space-time evolution of sec- 
ondary electron cascades in diamond. These cascades follow the impact of a single primary electron at energies 
between 0.5 — 12 keV, representing the usual range for photoelectrons. The calculations describe the secondary 
ionizations caused by these electrons, the three-dimensional evolution of the electron cloud, and monitor the 
equivalent instantaneous temperature of the free-electron gas as the system cools during expansion. The dis- 
sipation of the impact energy proceeds predominantly through the production of secondary electrons whose 
energies are comparable to the binding energies of the valence (40 — 50 eV) and the core electrons (300 eV) in 
accordance with experiments and the models of interactions. The electron cloud generated by a 12 ke V electron 
is strongly anisotropic in the early phases of the cascade (i < 1 fs). At later times, the sample is dominated by 
low energy electrons, and these are scattered more isotropically by atoms in the sample. The results show that 
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the emission of secondary electrons approaches saturation within about 100 fs, following the primary impact. 
At an impact energy of 12 keV, the total number of electrons liberated in the sample is < 400 at 1000 fs. The 
results provide an understanding of ionizations by photoelectrons, and extend earlier models on low-energy 
electron cascades (E = 0.25 keV, [Jl|,g]) to the higher energy regime of the photoelectrons. 



In atomic or molecular samples exposed to X-ray radiation damage occurs. In light 
elements it proceeds mainly via the photoelectric effect. Photoelectrons and Auger electrons 

are then emitted. They propagate through the sample, and cause further damage by 
excitations of secondary electrons. Photoelectrons released by X-rays of ~ 1 A wavelength 
(~ 12 keV) are fast, v ~ 660 A/fs, and they can escape from small samples early in an 
exposure. In contrast, Auger electrons are slow (v ~ 95 A/fs in carbon), so they remain 
longer in a sample, and it is likely that they will thermalize there. An analysis of electron 
cascades initiated by Auger electrons in diamond is described in |J2p. 

A detailed description of electron cascades initiated by an electron impact of energy be- 
tween ~ 0.5 — 12keV is needed for a better understanding of radiation damage in larger sam- 
ples as secondary ionization caused by propagating photoelectrons is significant there. Such 
processes need to be investigated for planned experiments with free-electron-lasers (FEL). 
Atomic and molecular clusters irradiated with VUV radiation have already been studied ex- 
perimentally 0]. 

Energetic photoelectrons (E > 0.5 keV in carbon) propagate almost freely through the 
medium. They interact with single atoms, and this interaction is well described by the Born 
approximation Here, we use the formalism of the Lindhard dielectric function [|7|] 

together with two optical models (the TPP-2 and Ashley's models p-|T3[]) to describe the 
inelastic interactions of electrons with atoms. This approach takes into account the valence 
and core ionizations of an atom, following the inelastic scattering of a free electron. This 
approximation works well for energies up to about 10 keV, and in this energy regime core 
ionization are responsible for not more than about 10 % of the total number of ionizations in 
solids [0. 

In this paper elastic scattering is treated in the muffin-tin potential approximation ||lj,]2|,]5|, 
^p. In previous studies on low-energy electrons (with energies up to, E = 0.4 keV), we used 
programs from the Barbieri/Van Hove Phase Shift package JT3|] . For electrons of energies, 
E = 0.4 keV, and higher, we obtained elastic cross sections from the NIST database [[16|]. 

Fig. [j] shows the total elastic and inelastic cross sections obtained from the calculations. 
The results show that for energies higher than 1 keV, the elastic and inelastic cross sections 
are comparable, but for lower energies, 0.1 < E < 1 keV, the elastic cross section is twice 
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as large as the inelastic one. For very low energies, E < 0.1 keV, the inelastic cross sections 
drop rapidly as the energy decreases, and elastic interactions become predominant. This was 
true for both optical models (Ashley and TPP-2). 

Fig. |2| shows the energy loss function (ELF) and the cross sections for core ionization 
from the K shell of carbon as estimated from the Lindhard approximation, using the optical 
models above. The figures were compared to results from the relativistic binary- encounter- 
Bethe model (RBEB) [Jl7|] for the total core ionization by impact electrons. This model 
was developed by combining a modified form of the Mott cross section and the leading 
dipole part of the Bethe cross section [p"8|]. Recently, it was extended to incident electrons 
with relativistic energies [[F7[]. When the incoming electron is fast, we can use the Bethe- 
Fermi approximation [ [T9| , [20] ] . It replaces the electric field of the incoming electron by an 
electromagnetic pulse of the same, short duration. Impact ionization is then proportional 
to the dipole transition probability caused by the short, non-periodic electromagnetic pulse. 
Incidentally, that is the basis of the two optical models used here. The ionization probability 
is proportional to the overlap of the bound electron wave function with that of the secondary 
electron. Simple scaling considerations predict that if valence electrons are ionized, their 
kinetic energies should not be much larger than a few times the shell binding energy. Figure 
3 shows that the mean energy of the secondary electrons in diamond is about 40-50 eV (close 
to the energy of L shell electrons), reaching peak energies of about 300 eV (close to the 
energy of core electrons). The energy loss function in Fig. ^is the sum of the large valence 
and small core contributions which add to the valence contribution at energies, E ~ 0.3 
keV. For comparison, the binding energy of the K shell in carbon is, Eb ~ 285 eV. It is 
difficult to separate the contribution of the core excitations from the valence excitations on 
the base of the ELF alone. Therefore, we make here a rough estimate of the pure core 
contribution by subtracting the valence component from ELF. The valence component was 
extrapolated at the edge of the core peak, i. e. for energies larger than 277 eV (cf. Fig. ||). The 
cross section obtained from this core contribution underestimates the RBEB cross section at 
energies smaller than 1 — 2 keV. This is due to the fact that the optical approximation does 
not apply at the core threshold. At larger energies our estimate agrees well with the RBEB 
predictions. This discrepancy does not make our results inaccurate. 

Fig. ^] shows the energy loss, to, during a single inelastic scattering of an electron on an 
atom as a function of the electron energy. These are the values obtained at fixed integrated 

duo' \da(E) / W) / / du>'(da(E)/du') , and repre- 

- — ! min J ^min 

sent the integrated probabilities of inelastic electron-atom scattering with the energy loss 
that is less or equal to uj at the impact energy, E. The energy loss, cu, does not exceed 0.3 
keV, and its value becomes independent from the impact energy, E, at energies greater than 
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1 keV. This is expected as one-electron excitations are predominant in inelastic scatterings 
0§- 



Results 

We performed a large set of Monte Carlo simulations in diamond, showing the path of an 
impact electron and the secondary electron cascade triggered by the electron. Electron tra- 
jectories were simulated as described in [0] in such a way that no energy loss was permitted 
to the lattice. This approximation gave a better estimate of ionization rates in [[|] than the 
approximation where the energy loss to the lattice was allowed. The starting position of the 
impact electron at t — fs was at the origin x=(0,0,0) of an arbitrarily chosen coordinate 
system, and the velocity of the impact electron was along the Z-axis. Motions of holes and 
ionizations by the holes were neglected since they influence the dynamics of electrons at 
very low electron energies only. The space-time characteristics of secondary cascades were 
recorded as functions of the impact energy and time. 

The energy dependence of the elastic and inelastic cross sections had a strong influence 
on the dynamics of the electron cloud (cf. Fig. [I]). At longer times, when most of the elec- 
trons have already cooled to low energies, elastic (isotropic) scatterings dominate the sam- 
ple. Electrons then propagate randomly through the sample, and their distibution is isotropic. 
Earlier studies show similar behaviour for Auger electrons [Q] . 

Evolution of the cascade was analysed through (a) the number of secondary ioniza- 
tions, N ei (t), and (b) the equivalent instantaneous temperature of the free electron gas kT(t). 
Quantities (a), (b) were averaged over a number of cascades. Figure |] shows the results ob- 
tained with Ashley's and the TPP-2 optical models for impact energies of E = 0.5, 1.5, 5, 12 
keV. 0.5 keV corresponds to the lowest energy of an incident electron at which the Born 
approximation may be applied when calculating the electron-atom cross sections in carbon 

The number of ionization events within the first femtosecond is within 10-25 for electron 
impact energies ofE = 0.5—12 keV. For impacts at lower energies, E < 1 .5 keV, the number 
of ionizations was between 20 — 70, and saturated within 100 fs. At higher impact energies 
there was an increase in the number of the ionization events after t > 100 fs, however, the 
sample was then dominated by low energy electrons (E < 20 eV) (cf. Figs. ^|), and thus not 
many ionizations could be expected to occur. We checked that at 1000 fs, there was a total of 
about 400 ionizations induced by an electron impact of 12 keV (Ashley's model), compared 
to 370 ionizations at 100 fs. When using the TTP-2 model, the corresponding numbers were 
~ 360 both at 100 fs and 1000 fs at 12 keV impact energy. 
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The equivalent instantaneous temperature of the electron gas decreased as the cascade 
evolved. We used this temperature equivalent since the electron gas was far from thermal 
equilibrium. However, the equivalent temperature is still a quantity conserved in electron- 
electron collisions. It is worth noticing that after ~ 10 fs the temperature curves obtained at 
totally different primary energies were very similar. This indicated that the number of elec- 
trons was not much influenced by the energy of the primary electron but rather by secondary 
electrons of lower energies (< 60 eV) which dominated the sample after 10 fs. At 100 fs 
the temperature of the electron gas dropped to ~ 5 eV (Ashley's) or ~ 2.5 eV (the TPP-2 
model). 

Fig. |5] shows plots of (a) the average number of electrons released, N e f, (b) the equivalent 
instantaneous temperature, kT, of the electron gas as a function of the energy, E, of the 
primary electron. These curves describe results obtained at different times (t = 1, 10, 90 fs), 
and were based on Ashley's model and on the TPP-2 model with no energy transfer allowed 
to the lattice. The data represent primary energies of E = 0.25, 0.5, 1.5, 5, 12 keV. Results 
for impact energy E = 0.25 keV were taken from Ref. [0]. The results can be used for the 
interpolation of the number of ionizations, and the temperature of the electrons at energies 
ranging from E = 0.25 keV to E = 12 keV. 

Fig. |5] shows that the number of secondary electrons, N el , is approximatly proportional 
to the impact energy. The slight curvature of the dependence indicates that the system has 
not yet reached equilibrium at 90 fs after the primary impact (cf. Fig. 0). 

The energy distribution of secondary electrons. The positions and velocities of elec- 
trons recorded at times, t = 1, 10, 90 fs at energies ofE = 0.5, 1.5, 5, 12 keV were collected 
from all cascades, and put into one file. Using these data, histograms for the energy distribu- 
tions were obtained, N(E)/N, at these time points. Number, N(E) = £f£°i Ni(E), is the 
average number of electrons in a bin, (E, E + AE), averaged over a number of cascades. 
Correspondingly, Ni(E) was the number of electrons found in that bin for the ith cascade. 
These distributions were normalized to the total number of electrons, N = J2e N(E). Fig.|| 
shows the histograms at impact energies of E = 0.5 keV and E = 12 keV. As expected, the 
energy histograms show that the number of low-energy electrons increased with time. One 
may notice that the dissipation of the impact energy was fast. At 1 fs, most electrons had en- 
ergies lower than 300 eV in the cascade. This follows from the assumption that one-electron 
excitations dominate electron-atom interactions, and events in which energy transfer exceeds 
the threshold for core ionizations by a secondary electron are rare. At 10 fs most of the elec- 
trons had energy lower than 40 eV in all cascades. The overall low energy of the electrons 
at that time influence ionization rates and the expansion rate of the electron cloud, and these 
both slowed down after 10 fs. At 100 fs there were only low energy electrons (E < 20 eV) 



5 



present in the sample, and elastic scatterings dominated the dynamics of the electron cloud. 

Spatial distribution of secondary electrons. In order to describe the spatial distribution 
of the secondary electron cloud in cascades triggered by E = 0.5 and E = 12 keV electrons, 
results from 500 simulations were analysed at these energies. Most electrons lay within a 
sphere whose radius depends on the electron impact energy. At 1 fs, it was r doud ~ 50 A 
for 0.5 keV impact energy, and r doud ~ 600 A for 12 keV impact energy. At 100 fs the 
corresponding radii were between 150 A and 3100 A respectively. At this time, the center 
of mass of the cloud moved 23 — 1400 A away from the starting point in the direction of 
the primary impact at energies of 0.5 — 12 keV, respectively (Fig. [7]). The shifts in other 
directions were small (1 - 33 A). 

N 



The sphericity tensor, S 
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, calculated for the collective cloud 



/ 

was almost diagonal. Non-diagonal elements were ~ 10 1 — 10 4 times smaller than the 
diagonal ones. The sphericity tensor was diagonalized, and the sphericity parameter, S, was 
obtained from its eigenvalues: S = 3/2 (A 2 + A 3 ), where Ai > A 2 > A 3 . For electron impacts 
of 0.5 and 12 keV S falls within the interval of 0.91 - 0.96 at 10 fs and 0.91 - 0.98 at 90 
fs, suggesting that the spatial distribution of the collective cloud was practically isotropic at 
that stage. 

Strong differences at the spatial distribution of the electron cloud manifested at 1 fs. At 
low impact energies (0.5 and 1.5 keV), S was within the interval of 0.8 — 0.9, while at 5 keV 
S was 0.4—0.5, and at 12 keV, S was around 0.30—0.35. This shows a strong anisotropy with 
increasing impact energy. The distribution of the cloud was cylindrically symmetrical along 
the primary impact vector but it was elongated along the direction of the impact. This is due 
to the fact that at 1 fs, the primary electron was much faster than the secondary electrons. 
At later times the sample became dominated by low energy electrons, which were scattered 
more isotropically. This is reflected in the progressing isotropy in the distribution. 

The final cloud at higher impact energies is pear-shaped. Figure Ogives a quantitative de- 
scription of the development of such clouds, using parameters shown in Figure |8|. Parameters 
z + , z_ and r represent root mean square values from 500 simulations, 



Z+ = \J {(z - Z C m) 2 ), Z > Z C M 
Z- = \J{{Z - Zom) 2 ), Z < Z C M 



r = \J((x- xcm) 2 + (y~ Vcm) 2 }, (1) 

and were calculated in respect to the center of mass of the cloud. The results confirm that 
the cloud is elongated along the axis of the primary impact (Z) at early times, and becomes 
more spherical later. 
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Figure 1 : Elastic and inelastic total cross sections for diamond. Inelastic cross sections are obtained from 
the Lindhard approximation with core ionization taken into account using the TPP-2 optical model or Ashley's 
optical model. Elastic cross sections up to energies, E — 0.4 keV, were derived with the Barbieri/Van Hove 
Phase Shift package fll5[]. For larger energies, the elastic cross sections were taken from the NIST database. 

Conclusions 

Our results can be used to estimate damage by photo electrons (E = 0.5 — 12 keV) in dia- 
mond and other carbon-based covalent compounds. The Monte-Carlo code may be adopted 
to simulate multiionization phenomena in different systems, ranging from the explosion of 
atomic clusters to the formation of warm dense matter and plasmas. The model could also be 
used to estimate ionization rates and the spatio-temporal characteristics of secondary electron 
cascades in biological substances. 



Acknowledgments 



We are grateful to Gyula Faigel, Zoltan Jurek, Michel A. van Hove for discussions. B. 
Z. thanks Peter Druck and the Institute of Theoretical Physics II at the Ruhr University in 
Bochum for providing the access to the computational unit. This research was supported in 
part by the Polish Committee for Scientific Research with grant No. 2 P03B 051 19, the EU- 
BIOTECH Programme and the Swedish Research Councils. A. S. is grateful to STINT for a 
distinguished guest professorship. B. Z. was supported by the Wenner-Gren Foundations. 



7 



References 

[1] B. Ziaja, D. van der Spoel, A. Szoke, and J. Hajdu. Phys. Rev. B, 64:214104, 2001. 

[2] B. Ziaja, A. Szoke, D. van der Spoel, and J. Hajdu. Phys. Rev. B, 66:0241 16, 2002. 

[3] R. Neutze, R. Wouts, D. van der Spoel, E. Weckert, and J. Hajdu. Nature, 406:752-757, 
2000. 

[4] Thomas Moeller et al.' Project on "Optical Non-Linear Processes, Multi-Photon Ioni- 
sation, and Coulomb Explosion of Atoms and Clusters by Intense Vacuum Ultraviolet 
Radiation". FEE, HASYLAB, DESY, 2000-2002. 

[5] B.H. Bransden and C.J. Joachain. Physics of atoms and molecules. Longman, Essex, 
pages 505-513, 1998. 

[6] PG. Burke. Atomic, molecular and optical physics handbook. Editor G.W.F. Drake. 
Woodbury, New York American Institute of Physics, page 536, 1996. 

[7] J. Lindhard. K. Dan. Vidensk. Selsk. Mat. Fys. Medd., 28:1-57, 1954. 

[8] S. Tanuma, C.J. Powell, and D.R. Penn. Surf. Interf. Anal, 11:577, 1988. 

[9] S. Tanuma, C. J. Powell, and D. R. Penn. Surf. Interf. Anal, 17:911, 1991. 
[10] K. Yoshihara S. Tanuma, S. Ichimura. Appl Surf. Scien., 100/101:47, 1996. 
[11] J.C. Ashley. J. Elec. Spec. Rel Phenom., 46:199, 1988. 
[12] J.C. Ashley. J. Elec. Spec. Rel. Phenom., 50:323-334, 1990. 
[13] J.C. Ashley. J. Appl. Phys., 69:674, 1991. 

[14] P. Pianetta. Low-energy electron ranges in matter. X-ray data booklet, pages 3-5, 2001 . 

[15] A. Barbieri and MA. van Hove, private communication; Phase Shift package, 
littp : //electron .lbl . go v/leedpack/ . 

[16] A. Jablonski and C. Powell. NIST Standard Reference Database. 

[17] Y.-K. Kim, J.P Santos, and F. Parente. Phys. Rev. A, 62:052710, 2000. 

[18] Y.-K. Kim and M.E. Rudd. Phys. Rev. A, 50:3954, 1994. 

[19] E. Fermi. Z. Phys., 29:315, 1924. 



8 



[20] H. Bethe. Ann. Phys. (Leipzig), 5:325, 1930. 

[21] J.M. Fernandez- Varea et al. Nucl. Inst. Meth. Phys. Res. B, 108:35-50, 1996. 



9 



a) 



b) 



10' 



1 tr- 



io 



10" 



10' ! 



-2 



10 



-8 



10 



10 




10 100 1000 10000 100000 

E [eV] 



RBEB (core) ASHLEY 
TPP-2 



o 

tT 



0.003 



eg 

^ 0.002 



0.001 



100 




1000 



10000 



100000 



[eV] 



Figure 2: (a) Energy loss function of diamond, 7m(— 1/e), and (b) the total cross section for core ionization 
of diamond. Results from the Lindhard approximation, obtained with Ashley's and the TPP-2 models are 
compared to the prediction of RBEB model for core ionization from K shell in carbon Jl7|]. The binding 
energy of a K shell electron in carbon is, Eb ~ 285 eV. 
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Figure 3 : Energy loss, u), in a single inelastic scattering event during electron-atom interactions in diamond 
as a function of electron energy. Plots show results at fixed integrated probabilities, P(u)\E) = 0.5, 0.8, 0.95, 
obtained with the TPP-2 model. The probability P(uj\E) is the integrated probability that the energy loss in a 
scattering is less or equal ui. Results from Ashley's model were very similar (not shown) 
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Figure 4: (a) Average number of secondary electrons emitted, N e i, vs. time; (b) The equivalent instantaneous 
temperature kT of electron gas vs. time averaged over 500 cascades. Curves correspond to the results obtained 
at different electron impact energies E = 0.5 — 12 keV from Ashley's model and the TPP-2 model with no 
energy transfer allowed to the lattice. 
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Figure 5: (a) Average number of secondary electrons emitted, N e i, vs. energy, E; (b) The equivalent instan- 
taneous temperature kT of electron gas vs. energy, E averaged over 500 cascades. Curves correspond to the 
results obtained at different times t = 1, 10, 90 fs from Ashley's model and the TPP-2 model with no energy 
transfer allowed to the lattice. The data were collected at primary energies of E — 0.25, 0.5, 1.5, 5, 12 keV. 
The results at energy E — 0.25 keV were taken from |Q| 
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Figure 6: Energy distribution, N(E)/N, (fraction of electrons per bin) among electrons (histogram) at (a) 
t = 1 fs; (b) t = 10 fs; and (c) t = 90 fs. Histograms correspond to results obtained at electron impact of 
E = 0.5 keV (left) and E = 12 keV (right) from Ashley's model (solid line) and the TPP-2 model (dotted line) 
when no energy transfer to the lattice is allowed. 
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Figure 7: Shifts of the center of mass of the cloud in the direction of the primary impact (Z-axis) vs. energy, 
E, in respect to the point of emission of the primary impact. Data on the electron cloud were collected from 
500 cascades. Curves correspond to the results obtained at different times t = 1, 10, 90 fs from Ashley's model 
and the TPP-2 model with no energy transfer allowed to the lattice at primary energies of E — 0.5, 1.5, 5, 12 
keV. 
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Figure 8: Schematic plot of the electron cloud. 
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Figure 9: Parameters z+, Z- and r describing the spatial structure of the electron cloud. The data from 
500 cascades were collected at times: (a) t = 1 fs, (b) t = 10 fs, (c) t = 90 fs at primary energies of 
E = 0.5, 1.5, 5, 12 keV from Ashley's model and the TPP-2 model with no energy transfer allowed to the 
lattice. 
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